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■ Abstract: 

^ . Quantum size effects for an exciton attached to a spherical quantum dot are 

^ i calculated by a variational approach. The band line-ups are assumed to be type- 

O I II with finite offsets Ve and Vh- The dependence of the exciton binding energy 

' upon the dot radius R and the offsets is studied for different sets of electron and 

^ hole effective masses. 
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I. INTRODUCTION 



With the development of new techniques for fabrication of semiconductor het- 
erostructures, quantum-size effects in these low-dimensional structures (quantum 
wells [|I], |, |, g |, |, 0], wires g and dots|, have been studied extensively [|n[ . 
The quantum state of an Wannier exciton is one main subject in this respect, 
since, by the spatial confinement in these microstructures, it is qualitatively dif- 
ferent from the exciton state in bulk materials. 

Several theoretical studies, with various degrees of sophistication, exist for ex- 
citons in type-I heterostructures for quantum wells, quantum wires and quantum 
dots. Excitons in type-II heterostructures have also been studied|P, |10|, p!3[] , 
but to a much lesser extent. A fundamental feature of type-II heteroj unctions 
is the spatial separation of electron and hole, which leads to longer radiative 
lifetimes, lower exciton binding energy, and unusual dynamic and recombina- 
tion properties of charge carriers as compared with type I. For an overview, 
including potential applications, we refer to a recent article |T^, from which we 
quote: "Though some III-V, IV- VI and II- VI semiconductor materials can form 
type II junctions (AlInAs/InP, InAsSb/InSb, InAs/GaSb, GalnAsSb/GaSb, In- 
GaAs/GaAsSb, Si/Ge, ZnTe/ZnSe, etc) the intriguing properties of these re- 
markable structures are still poorly understood." 

The case of excitons in type-II quantum dots has only recently been considered 
in two model calculations. Rorison[lT5| has used a simple separable wave function 
in a variational calculation, with parameters appropriate for GaAs/AlAs and 
InAs/GaSb dots. We[|16[ have, on the other hand, used a more sophisticated 
variational wave function, and also presented analytical considerations of limiting 
cases, with the aim to obtain insight in how the exciton binding energy and the 
electron-hole correlations depend upon the effective masses and the dot radius. 
In Ref. |]T6[ the offsets were assumed to be infinite, corresponding to a complete 
spatial separation of the electron and hole. It was shown that two different 
regimes exist: For dot radii R much smaller than the effective bulk Bohr radius 
the electron and hole are essentially uncorrelated, while for R much larger than 
the Bohr radius the electron and hole are strongly spatially correlated, residing 
near the dot boundary just opposite each other. 

This infinite barrier model is artificial for small dot radii R, since in realistic 
situations the confined particle then tends to leak out in the barrier material. 
Moreover, the importance of the incomplete confinement in optical experiments is 
very clear, since the oscillator strength for excitonic transitions is proportional to 
the square of the electron-hole wave-function overlap. Thus it will be worthwhile 
to make clear the dependence of the exciton binding energy, and of the wave 
function overlap, on the magnitudes of the offsets. For estimates of band offsets 
and effective masses for several III-V heterostructures see Refs. IT^ |TB|. 

In this work, we report the result of a variational calculation of the binding 
energy for the electron-hole system with finite offsets. Since it is always useful to 
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have results for limiting cases to compare with, we compute first binding energies 
for the case in which one particle is confined, and the other completely free, i. 
e. Vh = 0. Since zero offset for one of the particles is a situation intermediate 
between a type-I and a type-II heterostructure, we denote this special case as 
being of type I^. A limiting type-l| situation occurs when the confined particle 
is completely confined within the dot {Ve = oo). In Sec. Ill the dependence of 
the exciton binding energy upon the dot radius, the offsets and the effective-mass 
values is studied, while Sec. IV contains results for the electron-hole overlap in 
the wave function, an important factor for the magnitude of oscillator strengths. 
We summarize our findings in Sec. V. 

II. MODEL AND METHOD 

In the effective-mass treatment the Hamiltonian for the electron-hole pair 
which forms the exciton is given by 

H= I Ph 

where 

Veir) 

and 

Vhir) 

Type-II situations correspond to Ve and Vh having the same sign. For definite- 
ness we assume both positive so that the electron is the confined particle. The 
alternative configuration presents of course the same computational problem. 

For simplicity we have assumed the same electron (and hole) effective mass 
in the dot and the barrier material, and that the dielectric constants of the two 
media can be accounted for by a single average value e. With degenerate valence 
band, we let rrih be either the heavy-hole or the light-hole effective mass, thus 
neglecting the complications due to the off-diagonal terms in the Kohn-Luttinger 
Hamiltonian ||19|| . 

The main task is now to determine the ground-state energy Eq of the Hamil- 
tonian ( |I|). The binding energy Ei, of the exciton is then the energy required to 
remove the hole, given by the difference between Eq and the confinement energy 
Ec, the ground-state energy of the electron in the spherical dot: 

Eb = E,- Eo, (4) 
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for r < i? 
Vp for r > R 



(2) 



Vh for r < R 
for r > R. 



(3) 
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As units of energy and length we use an effective Rydberg energy 
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2r V47re 
and an effective Bohr radius 



ah = (6) 



both in terms of the hole mass. In terms of these quantities we denote 



as our dimensionless energy, and 



(7) 

^Ry 



R = ^ (8) 



as the dimensionless radius. 

The ground-state energy of the Hamiltonian ( [|) is determined variationally, 
using the following nonorthogonal basis set of functions: 



V'(re, r,) = ^ ^ Cm e-'^'^^ e"'^^''' g-^l^-^'^l (9) 

k=l 1=1 

The variational parameters are 1], 7, j3, and the expansion coefficients Cm, al- 
together 3 + NeNh parameters. Optimalization with respect to the expansion 
coefficients Cm is a generalized eigenvalue problem, consisting in diagonalization 
of a matrix of size N^Nh x NeNh- All matrix elements can be evaluated analyti- 
cally in terms of error functions. The only numerical problems that need special 
attention arise because the analytic expressions contain terms that individually 
diverge when + 1i — 2/3 is close to zero, or because exponential terms may 
become inadmissably large for small 77s. Finally, the minimalization with respect 
to 7], 7 and (3 is demanding due to the existence of many local minima. 

Good accuracy is usually achieved with a moderate size of the basis set. This is 
illustrated by Table 1, which shows, by means of a special example, how the result 
for the ground-state energy of the exciton depends upon the size of the basis set. 
Apparently a basis set of 9 functions (A^e = Nh = 3) suffices to give the energy 
within 1% accuracy. The exceptional cases occur for very large offsets, because 
one cannot achieve an almost vanishing wave function in the nearly inaccessible 
region with a small number of basis functions, and for large dot radius R when 
the exciton is located near the boundary. In these special cases we make use of 
alternative functional forms, as discussed below. 

One may of course also check the quality of the Gaussian basis set for the 
confined particle by comparing the one-electron energies computed variationally 
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with the known exact ground state. With = 5, the accuracy is in general on 
the 1% level, or better. Although the confinement energy in equation ( ^ may 
be obtained exactly, we compute it variationally for better procedural consistency. 



III. EXCITON BINDING ENERGY 

We first discuss the case of one particle confined to the dot material, with 
the other free to move. We denote this situation, intermediate between a type-I 
and a type-II heterostructure, as type Two examples of heterostructures with 



very small valence-band offsets are InAs/AlSb and ZnSe/Zno.79Mno.2iSe [^. 

A negligible hole offset, and a completely confined electron constitutes the 
extreme case: 

A. K = oo, V/, = 

For this extreme case we use, for the reason mentioned above, an alternative 
variational function 

sinfvrrp/i?) ^ i , s. 

re 

for Te < -R, zero otherwise. The numerical results for the exciton binding energy 
E\, are shown in figure 0. Fig. |^ exhibits the size effect of the three different 
contributions to the total energy, viz. the kinetic energies of the electron and 
hole and the Coulomb interaction energy. We see that for R smaller than about 
3 the electron kinetic energy is essentially equal to the ground-state energy in the 
dot, and dominates the other energy contributions. For large R the three energy 
contributions approach the ratio l:l:-4, characteristic of the two-particle problem 
in bulk. 

It is possible to understand the large- and small- i? behavior of the binding 
energy in figure [^. For large radii the binding energy must approach the bulk 
binding energy with a reduced mass: 

E,{R ^ oo) = . (11) 

rrie rrih 

For small enough dot radius, or small electron to hole mass ratio {R^rrif^/mh « 
1), the electron will reside in its ground state 'i/'o('"e) = sin(7rre/-R)/(reA/27r_R). 
Then the hole will see a charge distribution p(r) = — e^/'Q(r), corresponding to a 
potential 

V(r) = 

J 47re r — rJ 
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— -r — for r > i? 

-4f?7iil-2f?sin^ + 2 / sin2(7rx)f] (12) 
for r < R 

The radial Schrodinger equation with this potential gives a ground state that is 
in good numerical agreement with the rUe/ rrih = curve of figure |l|. For a finite 
effective-mass ratio this treatment can legitimately be used for small R only, in 
which case the potential ( |l^) is merely a small perturbation of the Coulomb 
potential. Using the Coulomb potential as the unperturbed potential and the 
difference V{r) + e^/Aner as the perturbing potential first-order perturbation 
theory yields 

E, = l-{l-i,)R' + OiR'). (13) 

For equal effective masses, it can be checked numerically that this is a reasonable 
approximation for R < 0.5. This range of dot radii is so small that the parabolic 
top is not visible in figure |l|. 

The spatial correlation function 

^ _ ((re - (re))(rfe - (r/,))) ^ (rer^) 
•\/((re-(re))^)((re-(re))2) ^{ri){rl) 

gives a quantitative measure of correlations. Figure |^ shows the results. For 
small me, as well as for a small dot radius, the particles are weakly correlated, 
since the electron is essentially in its ground state. And, as could be expected, 
the correlations are larger in the present case than for an infinite hole offset. 

For larger electron mass the correlations increase. The limiting case of a very 
large mass ratio {rrie » mh) can easily be treated in the Born-Oppenheimer 
approximation, which yields the limiting behavior p2[ 



R 

C{me/mh -» oo) -> (15) 
v/i?2 + 9/(l-|vr-2) 

The correlation function for the largest mass ratio in figure 
B. Ve finite, = 

For finite electron offset, still in a type-l| situation, we obtain more accurate 
results by means of the variational function ( |^). The results for equal effective 
masses are given in figure have a minimum size in order to be able to bind the 
electron after the break-up of the exciton. The presence of the hole, however, 
which in type-l| situations merely is a satellite to the electron, makes it possible 
to have the exciton attached to smaller dots, dots for which the offset potential 
is insufficient to bind the single electron ||23[|. 



The figure shows that the binding is weaker when the dot is large, and when 
the offset is low. All graphs show that for a given offset there exists a radius 
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Rm for which the binding energy is maximal, and Rm increases with decreasing 
Ve- This trend can be understood by envisaging the confined particle to create 
a charge distribution, with which the freely moving particle then interacts. The 
ground state in a smeared charge distribution deviates from the Coulomb value 
increasingly more the more extended the charge distribution is since the overlap 
decreases. Our interpretation is that the maximum binding energy of the two- 
particle system corresponds to the minimum extension of the wave function of 
the confined particle. Since the one-particle ground state in a spherical dot is 
a well-known textbook example, one can calculate its width (r^)2 exactly. The 
width is very large both when the dot is very large, and when the dot is so small 
that the potential barely binds the electron. Thus a minimum width of the one- 
particle ground state must exist for a definite dot radius r^, and a straightforward 



calculation [22 1 yields the exact behavior 




(16) 

with a numerical constant c = 2.67. For the offsets Ve = 1,5 and 25 in figure ^the 
binding energy is maximal at dot sizes Rm = 2.35, 1.11 and 0.53, respectively. 

This corresponds to the values 2.35, 2.48 and 2.65, respectively, for Rm\f^e, 
showing that the interpretation of the maxima makes sense, for large offsets even 
quantitatively. 

The binding energy also depends upon the effective masses. As shown in 
figure ^, where Ve = 25, the binding energy decreases when the electron effective 
mass becomes smaller. The position of the maximum of Ei, is seen to increase 
with decreasing effective mass ratio nie/mh in close accordance with the square- 
root dependence of Eq. ( |l^). Numerically, the maxima in figure § occur at 
dimensionless radii 1.63, 1.11, 0.71, and 0.53 for the mass ratios nie/mh = 0.1, 
0.2, 0.5, and 1, respectively, while Eq. ( [1^) yields the values fm = 1.69, 1.19, 
0.76, and 0.53 for these mass ratios. 

Let us consider a specific case. As mentioned above InAs/AlSb is close to 
a type l| heterostructure. We quote from Ref. |2^: "Since the valence band 



offset of 0.11 eV between InAs and AlSb is so small, we may approximate it as 
zero, The electron, on the other hand, is well confined in the InAs material. 



since Ve = 1.37eV [|T^]. Using||2^ an average dielectric constant e ~ 13.6eo, an 
average light-hole effective mass rrih^ — 0.07mo, we find values of the dimension- 
less electron offset Ve to be more than 250, i. e. effectively infinite. Taking the 
average heavy-hole effective mass to be ruh^h — 0.5mo, we find Ve ~ 35, also very 
large. The relevant figure is thus Fig. 1. With the electron effective mass in InAs, 
rrie = 0.022mo, the mass ratios me/rrih are about 0.3 and 0.04, respectively. One 
thus have to interpolate between the mass ratios and 0.5 graphs in Fig. 1. (For 
the heavy-hole exciton this can only be used for dimensionless dot radii larger 
than about 2.3, according to the expression ( |T^). The lower curve in Fig. 5, 
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although corresponding to an offset 25 rather than 35, and a larger mass ratio 
me/rrih, gives also a rough idea of how the binding energy for the heavy-hole 
exciton depends on the dot radius.) 

C. Ve finite, Vh finite 

We finally investigate the effect of a finite hole offset. The binding energy as 
function of the quantum dot radius R is shown in figure Comparison with figure 
1^ shows that the binding energy is smaller when the hole offset is nonzero, as 
could be expected. We see, moreover, that offsets larger than about 15 yield 
essentially the same binding energy as infinite offsets. 

The figure corresponds to mg = m/j. The binding energy for equal electron and 
hole offsets is, however, very insensitive to the mass ratio me/rrih'- Decreasing the 
mass ratio from 1 to 0.5 lowers Eb at most by a few percent. This insensitivity 



was also found for infinite offsets [|T6|. However, for a large disparity between 
the electron and hole offsets, the binding energy is more sensitive to the mass 
parameters, as witnessed by figure ||. 

The figure also shows the interesting feature that with increasing values of V^, 
Ef, now decreases, while it increased'm the absense of the hole offset (Fig. |). How 
can this be explained? An increase of any of the two offsets is a positive definite 
perturbation on the Hamiltonian ( |I|), and leads therefore to an increase in its 
ground-state energy Eq, i.e. gives a negative contribution to the binding energy 
E^. Hence E^, would be lowered if the hole offset were increased for fixed An 
increase in the electron offset, however, has the additional effect of increasing the 
electron confinement energy E^. in equation ( |^), and the competition between 
the two effects determines the change in the binding energy Ef, = E^. — Eq. The 
variation of the confinement energy Ec is the dominating contribution when Vh = 
(Fig. ^), but cannot outweigh the combined influence on Eq from increasing 
both the hole and electron offsets by equal amounts (Fig. ||). 

For type-II heterostructures the offsets are equal if the bandgaps in the two 
materials are equal. For GaAs and InP the bandgaps are (at 300K) 1.42 eV and 



1.34 eV, respectively The dielectric constant in both materials is close to 



e = 12.7eo. The band offsets are small, estimated p5[| to be of order 0.2 eV and 



0.3 eV. Let us assume that this may be approximated by a common offset value 
of 0.25 eV, and let us assume further that the complications due to to strain can 
be neglected. With an average electron effective mass mg ^ O.OTmo, and average 
light-hole effective mass nihj — O.lOmo, and the average heavy- hole effective mass 
mh^h — 0.54mo, we obtain the very rough estimates 30 (light hole) and 5 (heavy 
hole) for the dimensionless offset values. Thus the exciton binding energy versus 
radius graph for offset value 5 in Fig. 6 is relevant for the heavy-hole exciton, and 
the graph corresponding to the light-hole exciton is squeezed between the two 
proximate graphs for offset values 15 and infinity. As noted above the results are 
insensitive to the value of the mass ratio. 
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IV. Oscillator strengths 



The hole-electron overlap is a decisive factor in determining the properties 
of the exciton in optical experiments. Since infinite offsets give zero overlap, 
the magnitudes of the offsets must clearly be extremely important in type-II 
situations. 



The oscillator strength of an optical transition is proportional to the factor p6 

2 

iP{t, r) d^r 



fo 



(17) 



Figure ^ shows how the overlap factor fo increases dramatically with decreasing 
offsets. This is expected, of course, since the particles are able to penetrate the 
dot boundary when the offsets are lowered. In figure || we illustrate a typical 
situation, for dimensionless offsets equal to unity, and dimensionless dot radius 
R = 8. The marginal electron and hole distributions, obtained from the two- 
particle wave function, are clearly centered on each side of, and away from, the 
dot boundary. The tails of these marginal distributions are small at the boundary, 
while the overlap, which corresponds to = rh, is maximal at, or very close to, 
the dot boundary. 

How the overlap factor depends upon the size of the quantum dot is more 
interesting. The electron wave function is squeezed out of the dot volume when 
the radius becomes small enough, just as it is for small offsets. Then the electron 
and hole are able to correlate, yielding a large overlap function. The decrease of 
the overlap function with increasing R, shown in figure ^ is in accordance with 
this. However, /q must finally increase with R as 



/o ^ c ■ AttR', R^ oo (18) 

where the constant c depends on the offsets and on the effective-mass ratio. The 
rationale behind ( |1^) is that for very large radii the exciton will, on the scale 
of the dot radius, sit very close to the boundary, so that the determination of 
the wave function is essentially determined through the solution of a plane-wall 
problem. The plane-wall problem determines the constant c, and the additional 
factor of AtcR^ in ( 0) comes from normalization and integration over the dot 
surface. For the exciton near a plane wall at z = we use the simple variational 
function '0(re,r/i) = {Ce - Ze){Ch + Zh)e~^\''''~'"'^ for Ze < Ce, Zh > -Ch, and zero 
otherwise. For the parameters of Fig. |^ we find for the variational parameters 
Ce = Ch = 0.85 and /5 = 0.15, which implies c = 1.62 ■ 10"''. The rise of fo 
according to ( |18[) is outside the range of radii shown in the figure. 

V. Summary 
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The present work is a model calculation in which we have, in the effective- mass 
treatment, performed a variational calculation of cxcitons in spherical quantum 
dots. We have in particular studied how the exciton binding energy Eh depends 
on the dot radius R, the effective masses rUe and ruh, and the offsets. The nature 
of the binding energy maximum, as function of the dot radius, is clarified. 
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Figure 1: The dimensionless exciton binding energy Eb for the offsets Vg = oo, 
Vh = 0, as function of the dimensionless dot radius R = R/ah- The values of the 
mass ratio rrie/mh are shown on the right-hand side. 

Figure 2: The electron kinetic energy (e), the hole kinetic energy (h), and the 
absolute value of the Coulomb energy (c), as functions of the dimensionless quan- 
tum dot radius R. These three contributions to the total energy are measured in 
the effective Rydberg (5). Here mg = m/j, Vg = oo, and Vh = 0. The dashed line 
corresponds to the ground state of the single electron. 

Figure 3: The spatial correlation function (30) as a function of the dimensionless 
radius R = R/ah of the quantum dot, for the case Ve = oo. The fully drawn 
curves correspond to Vh — 0, the dashed curves to = oo [16]. The two upper 
graphs correspond to mg = lOOm^i, the two lower ones to mg = m^. 

Figure 4: The dimensionless exciton binding energy Ef, for V/j = and several 
values of T4, as function of the dimensionless quantum dot radius R. The values 
of Vg are given on the graphs. Here mg = m^. In the variational calculation 
Ng = 5 and Nh — 2 have been used. The dashed line corresponds to T4 = oo, 
calculated less accurately with the simpler function (10). 



Figure 5: The dimensionless exciton binding energy Eh for Vh — and Ve = 25, 
as function of the dimensionless quantum dot radius R, for several values of the 
effective mass ratio m^/nih (shown on the right-hand side). In the variational 
calculation = 5 and Nh = 2 have been used. 

Figure 6: The dimensionless exciton binding energy E^ for equal effective masses 
and equal offsets, as function of the dimensionless quantum dot radius R. The 
offset values are shown on the graphs. In the variational calculation iVg = 3 and 
Nh = 5 have been used. Results corresponding to = V/j = oo have been taken 
from Ref. [16] and are shown as a dashed graph. 

Figure 7: The overlap factor (17) as function of the offset. Here mg = nih, 
Ve = Vh = V, and R = 3. In the variational calculation Ne — 3 and Nh = 5 have 
been used. 
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Figure 8: The radial probability densities 47rrfV'(rj)^, i = e,h for the electron 
position (e) and the hole position (h), as functions of the dimensionless distance 
from the dot center. Here R — 8 (indicated by an arrow), Wg — ''^h ^-nd = 
Vh — 1. The dashed line shows 10^47rr^'?/'(i"i ^Y- (The small wiggle near f = 3 is 
presumably due to numerical inaccuracies.) In the variational calculation Ne — 3 
and Nh = 5 have been used. 



Figure 9: The overlap factor (17) as function of the dimensionless dot radius 
R. Here rUe = mh and Vg = Vh = I. In the variational calculation Ne — S and 
Nh = 5 have been used. The overlap factor will eventually increase, for very large 
R according to (18). 



THE FIGURES CAN BE OBTAINED FROM THE AUTHORS 



13 



TABLE 1. Exciton energy Eq for different basis 



set sizes. Here R = 2, rrie = mh, and = Vh = 5. 



Ne\Nu 


1 


2 


3 


4 


5 


1 


1.505 


1.406 


1.381 


1.377 


1.376 


2 


1.467 


1.368 


1.343 


1.339 


1.338 


3 


1.460 


1.361 


1.336 


1.332 


1.331 


4 


1.457 


1.359 


1.333 


1.329 




5 


1.456 


1.357 


1.331 
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Addendum 

1. Derivation of Eq. (15): Correlations for large mass ratio 

We seek the correlation function (14) 

C = (19) 

\rl){rl) 



for the very special case of Ve = cxd, V/i = and rrie >> ruh- 

According to standard Born-Oppenheimer procedure one first solves the ground- 
state Schrodinger equation for the light particle (h), for fixed position ff. of the 
heavy particle. This is a pure Coulombic problem yielding 

with r = r/i — r^.. The eigenvalue for this problem is — e^/Svrea/i. 

The second step in the Born-Oppenheimer procedure is to solve the Schrodinger 
equation for the heavy particle, with the eigenvalue for the light-particle problem 
serving as an additional potential. Since this eigenvalue is merely a constant, the 
confining potential determines the wave function completely. The ground state 
is 

The correlation function (1) is now easily evaluated using r instead of f^. With 
rh = 'f^e + 'f^^^ have 



{rl) 



since (r) = by symmetry. The integrals involved are simple: 



(20) 



2\ _ C r^e-'^'l"^ dr _ ^ 



and 

A , _ - , - 

J^sm^TirjR) dre Vs 2n^ , 

Inserting these averages into ( ^0]), and using R = R ah v^e obtain 

c ^ 



/ 2\ _ I^rlsm^inre/R) dr^ _ ^ (]_ 

y el rR ■ 2/ / T-.\ 7 ~ o 



i?2 + 9/(1 - 3/27r2) 
15 



which is equation (15) in the text. 



2. Derivation of Eq.(16): The minimum width of the one-particle 
ground state 

The unnormahzed ground-state wave function ilj[r) for one particle of mass 
m in a dot potential 



is given by 

Here 
and 



sinA;e/^(^-^)/^ r > R. ^^^> 



k = Rj2mE/h^, (23) 



(3 = Rj2m{Ve - E)/Ti^ = J 2mR'^Veh~'^ - (24) 



with E the ground-state energy. 

The wave function ( E2|) is continuous. Continuity of %l)'[r) requires 



(3 = -kcotk. (25) 
The width < >^ of the one-particle ground state is thus given by 

R oo 

/ srn^ikr / R)dr + e^^ sin^ k J r'^e~'^^'^/^dr 



R oo 

/ sin^ {kr/R)dr + e^^ sin^ k J e-'^^'^/^dr 



R 

1 _ sin(2fc) _ cos(2fc) , sin(2fc) , sin^ fc ( I \ J_ \ ]_ 

_ p2 6 ik 4fc2 8A:3 /3 V 2 2/3 4/3^ y , - 

~ 1 _ sin(2fc) , sinifc ^^^> 

2 4A; ~'~ 2/3 

Using ( |2^) and ( ^) we have 

= ^ . (27) 
2'mVe sin A; 

Thus we may express the width entirely in term of the variable k: 

2 ,2 1 _ sin{2k) _ cos(2fc) , sin(2fc) _ sin^ k (l _ 1 , 1 \ 

6 4fc 4fc2 8fe^ fccotfc V2 2k cot k ~^ Ak^ cot^ k J 



< > 



2mVe sin^ i - ^HiM _ ^iiiiL 

2 4fc 2k cot k 

{2i 
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Thus 



(29) 
(after 



\ 2mVe ■ 

where c is a number, is the minimum value of the function of A; in ( 
the multiphcation sign). Numerically we find c = 2.67. 

Finally, by introducing an effective Bohr radius as unit of length and an 
effective Rydberg as unit of energy, both in terms of the hole mass, equation 
791) takes the dimensionless form 



(30) 




This is equation (16) in the manuscript. 
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